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We study the effect on the primordial cosmological perturbations of a sharp transition from infla- 
tionary to a radiation and matter dominated epoch respectively. We assume that the perturbations 
are generated by the vacuum fluctuations of a scalar field slowly rolling down its potential, and that 
the transition into the subsequent epoch takes place much faster than a Hubble time. The behaviour 
of the superhorizon perturbations corresponding to cosmological scales in this case is well known. 
However, it is not clear how perturbations on scales of and smaller than the Hubble horizon scale 
at the end of inflation may evolve through such a transition. We derive the evolution equation for 
the gravitational potential ^ , which allows us to study the evolution of the perturbations on all 
scales under these circumstances. We show that for a certain range of scales inside the horizon at 
the end of inflation, the amplitude of the perturbations are enhanced relative to the superhorizon 
scales. This enhancement may lead to the overproduction of Primordial Black Holes (PBHs), and 
-) \ therefore constrain the dynamics of the transitions that take place at the end of inflation. 

£ 

H- 8 



o 



43 : 

Or 



(N 
> 
OS 



PACS numbers: 04.70.-s, 98.80.-k 



I. INTRODUCTION 



Inflation has become the main paradigm to understand the presence of small inhomogeneities during the early 
universe. The observed Cosmic Microwave Background (CMB) and Large Scale Structure (LSS) spectrum of pri- 
mordial perturbations, reveals that simple inflationary models are consistent with the analysis of the current data. 
Nevertheless, CMB and LSS observations only probe primordial perturbations on scales that leave the horizon long 
before the end of inflation, and possible cosmological implications of the perturbations on subhorizon scales have not 
been studied much. 

In this paper, we investigate the behavior of curvature perturbations on scales of and smaller than the horizon size 
at the end of inflation. We assume standard slow-roll inflation and that the perturbations are generated from the 
t— > . vacuum fluctuations of the inflaton field. It is expected that the universe enters either a stage of radiation-domination 
or matter-domination right after inflation, depending on the efficiency of preheating or reheating, and the transition 
takes place within a timescale shorter than the Hubble time. We therefore model this transition as a change of the 
equation of state from w « — 1 to w = 1/3 or to w = within a time scale At <C H^ 1 where H e is the Hubble 
parameter at the end of inflation. 

On large scales where the wavelengths are much larger than At, we may regard the transition as instantaneous, 
occuring on a given spacelike hypersurface. Then we can match the perturbations between inflation and the succeeding 
. epoch by requiring that the intrinsic and extrinsic curvatures are continuous on the hypersurface at which we do the 
CsJ 1 matching It is natural to assume that the end of inflation occurs when the inflaton field reaches a critical 

value, 4> c (x,t) = constant. This implies that the transition occurs on a comoving surface on which tf> is spatially 
« ■ homogeneous. Then one can show that we have the following junction conditions, 

On i 

O ■ [K k ]± = 0, (1) 

=0, (2) 

where 1Z is the curvature perturbation on comoving hypersurfaces, and $ is the generalized gravitational potential. The 
^3 ■ junction conditions give the initial conditions for the following radiation- or matter-dominated stage. On superhorizon 
scales, the potential $ then reaches a constant value, ~ 1Z, a few Hubble times after the end of inflation 

What we are interested in here is perturbations on scales of or smaller than the horizon scale. On subhorizons scales, 
but still larger than the scale At, the resulting perturbations are the same as in the case of superhorizon perturbations 
if the universe is matter-dominated after the transition. On the other hand, the result is not so simple if the universe 
is radiation-dominated after the transion. In Refs. (HQ, we investigated this case, and argued there that the relative 
large amplitude of the potential VE' can be achieved on subhorizon scales and may lead to the overproduction of PBHs 
at the end of inflation. The constraint on the amplitude of the curvature perturbation is then tighter than the current 
bound from the PBHs formed throughout the radiation epoch [5[ (for the constraints on the " standard" formation of 
PBHs see [HljJ). 
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In this paper we investigate the effect of the transition on curvature perturbations on all subhorizon scales, including 
those on sufficiently small scales on which the transition cannot be regarded instantaneous any more. The paper is 
organized as follows. In Sec. (jll]) we derive the governing equation for the evolution of the potential ^> for an arbitrary 
equation of state and velocity of sound. We review the generation of the perturbations during slow-roll inflation in 
Sec. piip. We study the evolution of VP through a transition into radiation and matter domination in Sees. (|IV|) 
and (O before concluding. In Appendices [X] and |Bj we present analytical approximations that we have employed to 
reproduce and understand certain features that we have found in our numerical estimations. 

Thoughout the paper, we use the Planck units 8nG = Mp 2 = 1. 



II. GAUGE INVARIANT COSMOLOGICAL PERTURBATIONS 



In this section, we introduce the governing equations for the scalar perturbations about the homogeneous and 
isotropic background that we will use through out this paper. 

The background space-time geometry for a spatially flat universe is given by the FLRW metric, 

ds 2 = a 2 (r) (-dr 2 + dx 2 + dy 2 + dz 2 ) , (3) 

where r is the conformal time defined in terms of the cosmic time t as dr = dt/a, and a(r) is the scale factor. Relative 
to the above coordinates the energy-momentum tensor of a perfect fluid at rest is 

n = -p, (4) 

Tj = PS], (5) 

where p and P are the energy density and pressure of the fluid. 

The governing equations for the system are determined by the Einstein Field equations and the energy momentum 
conservation. That is 

R^--g" v R = T^, (6) 
T" u ]V = 0. (7) 

In the present case they give the following set of equations: 

m 2 = pa 2 , (8) 
p + m (p + P) = , (9) 

where a prime ( ') denotes the differentiation with respect to r, and the conformal Hubble parameter % is defined 
by H — a' /a, hence is related to the physical Hubble parameter H by H = aH . A combination of the above two 
equations gives a second order differential equation for a(r), 

W = - l -(p + iP)a 2 . (10) 
6 

Ignoring the vector and tensor perturbations, the perturbed metric can be written to first order as 

ds 2 = a 2 (r) [- (1 + 2A) dr 2 - 2B^drdx l + [(1 + 2D) 5 %] + 2E tij ] dx l dx j ] , (11) 

where a comma denotes the partial differentiation with respect to the spatial coordinates, and A, B, D, and E are 
the functions that describe the scalar perturbations. The scalar part of the perturbations in the components of the 
energy momentum tensor for a perfect fluid are 

6T° = -Sp, (12) 
ST? = (p + P)(B, i + v, i ) , (13) 
ST] = 5P8), (14) 



where we have decomposed the fluid 3-velocity v l into a divergence free vector and the gradient of a scalar, v l , to 
extract the scalar part of the perturbations as usual. 
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It is convenient to express the governing equations in terms of quantities which are invariant under infinitesimal 
coordinate transformations. Bardeen [ItJ gives the following two gauge invariant quantities for the metric, 



$ = 



A+-[a(B- E')]' , 



* = D + U{B - E') 



and for the matter variables, 



6c 
v 



5 + 3^-^jH(v-B) , 
v-E' . 



(15) 
(16) 

(17) 
(18) 



The functions "J and $ represent the lapse function perturbation and the curvature perturbation in the longitudinal 
gauge respectively. For the matter variables 5c is the density contrast 6p/p on the comoving slice and v s the velocity 
of matter on the longitudinal gauge relative to the x l = constant observers [18j . 

For the matter perturbations of scalar type, there is a third gauge invariant quantity often called the entropy 
perturbation, which can be constructed in terms of the matter variables only: 



P' 

5P = 5P -Sp . 

P 



(19) 



The entropy perturbation measures the difference between uniform density and pressure hypersurfaces. For an equa- 
tion of state in which the pressure is a function of the energy density alone (the so-called barotropic equation of state) , 
the uniform density and pressure hypersurfaces coincide and the pressure perturbation is called adiabatic. For a more 
general equation of state this is no longer the case and we get a non-zero entropy perturbation. However, there is one 
important exception in which the perturbation is purely adiabatic although the equation of state is not barotropic. 
That is the case of a single real scalar field. In this case, P = K — V(<f>) and p = K + V(<fi), where K = 4> 2 /2. 
Then, it is known (or easily seen from the equations below) that one obtains a closed single second order differential 
equation for the scalar perturbation if one defines the sound velocity as the ratio between SP and Sp on the comoving 
hypersurfaces (denoted by the suffix C), 



SPc = c\Spc ■ 



(20) 



In the present case, c s = 1 since SPc — Spc — 6K. Note that, since Sp = c 2 Sp in any gauge for a barotropic fluid, the 
definition of the sound velocity c s by Eq. (|2U1) is valid both for a scalar field and for a barotropic fluid. 

Now we can express the perturbed equations for Eqs.© and ([7]) in terms of the gauge invariant variables above as 



V 2 * 
$ + * 

1 ( -v 

a 



pa 



0. 



--(p + P) a 2 (av) 



(P + P) 



5 Pi 



c ■ 



(21) 
(22) 
(23) 

(24) 



Then setting SPc = c\pSc in the last equation, we can reduce the equations above to the following second order 
differential equation for 



a 2 (p + P) 
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a 2 (p + P) 



+ 



, 2 2 

k c. 



\{P + P)a 2 



fk = 0, 



(25) 



where fk = a^k- The equation above determines the evolution of the potential ^ for a given equation of state and 
velocity of sound c s (t). In the following sections we use Eq. (|2"5"]) to calculate the evolution of the perturbations 
generated during slow-roll inflation through a rapid transition into radiation domination and matter domination 
epochs. 
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III. THE GENERATION OF PERTURBATIONS DURING INFLATION 



We assume the universe is dominated by a minimally coupled real scalar field (f> during inflation. As mentioned in 
the previous section, for this universe, the energy density and pressure of the homogeneous scalar field (f> = 4>(t) are 
given by 

p = K + V(<P), P = K-V(<P); K~ l -^ = ^ 2 . (26) 

The density and pressure perturbations are then related by 

SP = Sp- 2V (k S(j), (27) 

which leads to SPc = Spc on comoving hypersurfaces and hence c s = 1. During this stage, the evolution equation for 
fk, Eq. (|23)l . may be written as 

fk -2jfL+[k 2 + n' -H 2 ]f k =0. (28) 
Introducing a new variable u given by 

/ r 6" \ 2k 
u k = f k cxp \- J ^drj = 2k f k( j>'- 1 = -j* k , (29) 

the differential equation for / can be expressed in the form, 

< + W 2 (fc,r) u k = 0, (30) 

where W 2 (k, r) is given by 

W* (k, r) = k 2 c 2 s + m j ff a 2 ; rn 2 eff a 2 = W - U 2 - (£) ' + (^) ' . (31) 

The coefficient 2k in the definition of u k is chosen so that the function a q k becomes a properly normalized mode 
function for an effective scalar field when quantized (22[. To be more specific, the field Q = aq becomes a conformally 
coupled scalar field with its effective mass-squared ni 2 ff defined in the above equation. 

We can find approximate analytical solutions for u k as follows. We first introduce the following set of slow-roll 
parameters 

«-A~4(T>', (32) 



2H 2 H 2 \a 

VH = "Hi = ~W ' (33) 



Then we rewrite the function W (k, r) as 



W 2 (k, r) = k 2 ~ H 2 [2e H - m + r) 2 H - e Hm + ^ j . (34) 

It is well know that the differential equations for the evolution of cosmological perturbations have analytical solution 
if €h and t]h are constant [13, Hl[ . Typically the slow- roll parameters are small, at least for the range of field values 
at which perturbations on scales constrained by observation are generated. Their time derivatives then, which are 
given by 

^-i H = 2e H {e-H ~ m) > ( 35 ) 
ti 

JjT)H = VH ( + CH + VH J , (36) 
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are small, and e# and r\n can be approximately taken as constant, as least for a certain number of e-folds. In that 
case the conformal time is given by 



da 
~a T H 



1 
~a~H 



enda 
a 2 H 



1 



1 



aH 1 — en 



(37) 



where for the time being we have set the integrating constant to zero. Writing down explicitly the time dependence 
in Eq. (f30|) . it becomes a differential Bessel equation, 



k 2 - 



u k = . 



where the order v is approximately given by 



v ~ - + 2e H - r)H 



In the limit k 3> aH, we assume that the fluctuations of the field are in the Minkowski vacuum, 

1 



-ikr 



(2k) 



1/2 



(38) 



(39) 



(40) 



Then the solution for Uk can be expressed in terms of the Hankel function of the first kind H^(x), and apart from 
an irrelevant constant phase factor, ^k can be written as 



^ H *(-kT?l 2 HV{-kT) 



2fc 3 / 2 

On scales well outside the horizon, k -C aH, the amplitude of ^ takes the following asymptotic limit, 

2"- 1 r (v) 1/3 _„ 

— — -e H (-fcr) -y- . 



* fc (r) 



/2pr(i/2) 



(41) 



(42) 



Evaluating the expression above at a given time, for example at the end of inflation, it can be seen that the scale 
dependence on superhorizon scales is given by the deviation of v from the value 1 /2. Nevertheless, the scale dependence 
then is commonly given by the small variation of the amplitude, (H 2 j <p)\k=aH , when perturbations on different scales 
leave the horizon during inflation. In this case one can write 



(r) 



r M 



/ 2pr(i/2) 



e H (1 - e H f 2 - [ ^ 



H' 



(43) 



k=aH 



For scales that remain inside the horizon during the inflationary epoch, \Jj does not display significant scale depen- 
dence even if the slow-roll parameters are relatively large. For kr > 1, ^ converges very quickly to 



*fc(r) 



1 H 2 



/2fc3 



-e H e 



-ikr 



(44) 



and only perturbations on scales for which fcr ~ 1, may vary respect to the scale invariant case. 

A special case, in which the slow-roll parameters defined above are constant and we get an exact analytical solution 
for "ffc, is power law inflation. In power law inflation a oc t q , and the slow- roll parameters give 



1 



en = m = 



(45) 



Then the solution for is given by Eq. (14"TT) with 



(46) 
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For large q, when the slow-roll conditions hold strongly, ^ on superhorizon scales becomes 1 

1 1 fH 2 " 



(t) = 



0. 



(47) 



k=aH 



On cosmological scales, which leave the horizon long before the end of inflation, it is common to describe the small 
inhomogeneties observed in CMB temperature fluctuations by the curvature perturbation on comoving hypersurfaces, 
7Z. The curvature perturbation 1Z is related to the Bardeen potential VP by 



n 



2U- l W + * 
3 (1 + w) 



(48) 



where w = P/p and (1 + w) = 2ejj/3 during inflation. 

For power law inflation, the derivative of \t can be expressed as 

n-^' k (r) = ^f 2 ^H {-kr) 1/2 [(-kr) (1 - e H ) (-kr) - (1 + e H ) (-kr) 

Then we find the following exact expression for 7Z: 

Uk (T) = ~2fc&? (1 ~ 6h) ( " fcT)3/2 H & ( ~ fcT) • 



(49) 



(50) 



On superhorizon scales, taking the asymptotic limit of the Hankel function, (—kr) — > 0, the curvature perturbation 
can be expressed as 



t- 1 ! 2 r(i/ + i) 

%/2fc3 r(3/2) 



(,+f/ 2 r- 1/2 (^) , 

V 9 / k =aH 



in agreement with [20]. In the limit q — >• oo, we get 



(51) 



(52) 



k=aH 



which is the standard result for a slowly varying potential 22j. For cr not too small this result is modified by the 
factor in Eq. ([5"Tj) . 

In view of the expressions derived above for 1Z and ty, it is clear that for small values of the slow-roll parameters 
the amplitude of "J is much smaller than 7Z on superhorizon scales. On scales well inside the horizon, k > aH, the 
curvature perturbation 1Z increases with the comoving scale as k/aH, while ^ is scale independent so the amplitude 
of is further suppressed relative to 1Z. It is natural then to neglect \& respect to 1Z during inflation if the slow-roll 
parameters are not too close to 1. 

As we stressed above, Eqs. (|4Tj) and (|50| are exact solutions for the perturbed quantities 1Z and VP generated during 
power law inflation. They provide approximate solutions for slow-roll inflation if we consider a suitable power law for 
the scale factor of the homogeneous background during a sufficiently short lapse of time. The approximation will be 
valid then for certain number of e-folds, depending on the variation of the slow-roll parameters. Here we are interested 
in the behaviour of the perturbations at the end of inflation, and therefore the power law solutions at that time may 
differ considerably from those at the time of horizon crossing for scales which are well outside the horizon at the end 
of inflation. Nevertheless, since we focus on the scales smaller than the horizon, this problem does not affect our 
analysis. 

In the following sections, we give approximate analytical expressions for ^ for a transition from inflation to radiation 
and matter domination alongside of numerical calculations. For the numerical estimation it is convenient to normalize 
the solutions of ^> with respect to the value of 1Z at horizon crossing at the end of inflation. Thus we define this value, 
An, as 



A 



■R 



r(q) 



H 2 



(53) 



k=H c 



1 The relation in Eq. 1471 holds barring an exact de Sitter background, which corresponds to <fi = 0. Note that we are considering the 
quantity (H 2 /if) small to be consistent with observation. 



7 



where r(q) is just the factor in Eq. (I5T|) . That is 



' (9) = 2 "" 1/2 rW ((/ + 1/2r " 1/2 ' (54) 



IV. THE TRANSITION INTO RADIATION DOMINATION 



In this section, we study the behavior of the potential ^ during a sharp transition into radiation domination from 
a slow-roll inflationary epoch. The evolution of the potential ^ from Eq. (|25[) is determined by 



fk ~ 



(p+p) 



2U 



f'k + 



(r)-- 2 ( P - 



P)a 2 



fk = 0. 



(55) 



It is necessary to do a constant time shift of the conformal time to match the scale factor between inflation and the 
radiation dominated stage. If the transition into radiation domination occurs much more rapidly than a Hubble time, 
then the scale factor at the end of inflation is approximately given by a c ~ 1/H c t , and we have to replace in the 
formulae for the perturbations during inflation t by r — tr, where tr is the constant conformal time shift. Neglecting 
e H in Eq. fl37|), we get t r = 2r c . 

One can in principle solve Eq. (|55[) numerically for a given equation of state and velocity of sound. Here we are 
interested in the qualitative behaviour of $ for a sharp transition into radiation domination, so we assume that the 
equation of state changes abruptly to reach its radiation domination value in some small fraction of a Hubble time. 
We can parametrise then the equation of state as P = w(t)p, where w(t) is a step like function with a certain width, 
At = 5/H c . During inflation (1 + w) = 2ejy/3, which is small quantity under slow-roll conditions. For a rapid 
transition into radiation domination it increases very sharply until it reaches the value 4/3, corresponding to the 
radiation equation of state P = p/3. 

To solve Eq. (l55j) for / = , it is convenient again to express it in terms of u, = 2(k/a)(p + P)^ 1 ^ 2 .fk- The 
evolution equation for u takes the form, 



where W 2 (k, r) is 



and A(t) is defined by 



< + W 2 (fc, r) u k = , 



W 2 (k, t) = k 2 c 2 (r) - AttG (p +p)a 2 - \a 2 (t) + l -A! (r) , 



Mr) 



(P + P) 



During the transition we may write W 2 (k,r) as 

W 2 (k,r)~k 2 c 2 (r)-V cS (r) , 

where V c g is approximately given by 



(56) 



(57) 



(58) 



(59) 



V cB (r) 



f 



2(1 + ™)' 



-H 



l + w' 



(60) 



neglecting all the terms proportional to T-L 2 . The terms proportional to the derivatives of w are sharply peaked 
functions, with an approximate amplitude of order H 2 /S 2 during the transition when the equation of state varies. In 
Fig. [TJ we show the time dependence of V e s(r) for an equation of state such that during the transition (1 + w) it 
behaves like 



(! + «>)= 3 



1 + erf 



He (t - T c 



(61) 



where we have ignored the small value proportional to en it has during inflation. We plot Voff for S = 0.01 and two 
different values of q. The height and width of the peaks shown in the figure depend on S as well as which is 
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evident from the expression of V e g in Eq. ([50]) . As we shall see, the variation of V c g with e# is responsible for the 
peculiar scale dependence of the spectrum of on subhorizon scales. On superhorizon scales, one gets the usual red 
tilt on the spectrum for power law inflation. 

Once the equation of state reaches its radiation domination constant value, the derivatives of w vanish and W 2 (k, r) 
becomes 



W 2 (k,r) 



k c 



(62) 



which gives a Bessel differential equation of order v = 3/2 for u when c s is constant. Note that the solutions for half 
integer values of v are just spherical Bessel functions. 

For scales k ^> H/8 then, the gradient term in W 2 (k, r) dominates. Assuming that c 2 (t) drops down from 1 to 1/3 
during the transition, an approximate estimate for VE^ on these small scales, k 3> aH/S, throughout the transition 
and radiation domination epoch, is given by the WKB solution, 



*fe (r) 



1 



(.P + p) 



1/2 



(2fc) 3/2 yfiJj)M* 



-ik f c s {r)dr 



(63) 



where we have recovered the Plank mass, M p — (SuG)^ 1 / 2 , so that the dimensionless nature of •\/2fc 3 |^'/ c | is explicit. 
Hereafter, we call this non-dimensional amplitude as the amplitude of \Jj at a comoving scale k. The solution in 
this regime represents plane waves with an increasing amplitude during the transition that starts to decay once the 
radiation conditions ensue. The maximum amplitude of the evanescent wave depends only on the energy scale of 
inflation, ~ H c /M p . 

In the opposite regime, k <C H c /S, the variation of the equation of state from the slow-roll regime is responsible of 
the growth of the perturbed quantity We can neglect the terms proportional to fh in Eq. (|25D . and therefore it 
becomes 



fk ~ 



which gives the following conserved quantity, 



(p + p)' 
(P + P) 



f'k 



+ 2H 



(P + P) 



constant . 



(64) 



(65) 



Integrating the equation above, we find that \E' does not change substantially during the transition if the transition is 
much faster than a Hubble time. In fact, the relative variation of the potential \& during the transition is roughly 



A* 



_s_ 



(66) 



Therefore taking to be constant during the transition on these scales appears to be a good approximation in the 
limit 5 < eh . 

Neglecting the small variation of "J and the scale factor during the transition, Eq. (|6"5"|) implies that the curvature 
perturbation on comoving hypersurfaces, 1Z, remains constant during the transition. Thus for a rapid transition from 
slow-roll into radiation domination, there is a certain range of scales inside horizon for which the transition may be 
taken as instantaneous. Using the junction conditions in Eqs.((T]) and we get the result found in Q, 



* fc (r) = 2K k (r e ) 



[{x — x e ) cos [x — x e ) — (1 + xx e ) sin (x — x e )] 



(67) 



where x = c s kr and x c = c s fcr e , with the velocity of sound given by c s = l/-\/3- On superhorizon scales, x c — > and 
x <C 1, VP reaches the constant value, 



* fc ~ --TZ k (t c ) 



(68) 



whereas the modes inside the horizon at the end of inflation undergo damped oscillations. 

We can check the qualitative behaviour of the perturbed potential 5', outlined in the preceeding paragraphs, by 
solving numerically Eq. (I55p for a definite equation of state which behaves like a step function with a given width. 
Using an error function for the variation of (p + p) during the transition, with a typical width At ~ 8/H c , we have 
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found that the numerical solutions agree very well with our analytical expectations in both regimes. In Fig. [TJ we 
show the evolution of VP for two different values of q in power law inflation with comoving wavenumber k <C H c /S. We 
take values of 5 such that S < 1/q. We plot the modes with a = 0.01, 0.1, 1 and 5, where a = k/a e H e = k/H c - Thus 
the former two correspond to supcrhorizon perturbations at the end of inflation, while the latter two to subhorizon 
perturbations. The superhorizon modes reach the value tyj- = 27?.fc(r c )/3 a few Hubble times after the end of inflation. 
For large values of q the scale dependence of these modes is small. As q decreases the deviation from a scale invariant 
spectrum in the superhorizon modes becomes more prominent. On subhorizon scales, k > T-L c , ^ oscillates with 
decaying amplitude once the radiation dominated epoch begins. The relevant quantity to estimate the abundance of 
subhorizon PBHs formed at the end of inflation is the first maximum of the oscillations, ^M(k) [HIE]- \I/m(&) increases 
with k, until the comoving scale approaches k ~ H e /8. ^yiik) then decreases and converges to the approximate value 
given m Eq. ([63]), * M ~ H c /M p , as shown in Fig. H 

The numerical value of ^M(k) on subhorizon scales for two different values of q is shown in the left panel of Fig. [5] 
The value of r for the first maximum is approximately given by tm — t c + ir'H c /2kc s . It is convenient to express 
\1/m(&) in H c /M p units, because in the limit k 3> He/S it approximately drops down to that value. For a given 
transition time, At ~ S/H e , the amplitude of increases as q gets larger on an interval about the comoving scale 
k r~j Ti.e/5, where the resulting spectrum has a broad resonance. As we show in Appendix [B] the resonance can be 
interpreted as the effect of V e g in Eq. ([59")) acting as a potential barrier for the incoming wave function. The higher 
the potential barrier is the larger is the maximum amplitude of the resonance. The width of the resonance is related 
to the transition time only, and it is of the order H. e /8. As k becomes much larger than He/S, i/»m tends to H c /M p 
independently of the values of q and 5. 

In the right panel of Fig. [3] we show numerical values for the peak of the resonance, which here we denote it as 
^RES) as a function of q for different values of the transition time. The solid lines represent the fitting curve, 

V2tf 1-res (<7, 6) = y/q [A + B V5 + C5] , (69) 

where A = 2.47, B = —4.74 and C = —1.58 in H c /M p units. The estimated value of peak of the resonance above 
is in agreement with the result found in 4] , where the case of an instantaneous transition into radiation domination 
was analyzed. Calculating the maximum of Eq. d57J), we get 0] 

V2tf f RES ~ 2V3 (1 - e H ) . (70) 

This is the same value we get in Eq. (I69|) if we take the limit 5 — > and express in terms of A-n for en small, 
noting that during inflation, 



6 V2 \ m p 



(71) 



We can therefore have a rough idea of the shape of the wide resonance if q and 6 are given. The height of the 
resonance is of the order of the curvature perturbation at the end of inflation, An , while the width is approximately 
given by 1/d. It is interesting to note as well that while on superhorizon scales a smaller value of q tends to enhance 
the tilt of the spectrum, it suppresses the amplitude of the resonance and in general of all the modes well inside the 
horizon about the comoving scale k ~ 8/T-L e . 

V. THE TRANSITION INTO MATTER DOMINATION 

If an epoch of matter domination succeeds the inflationary epoch, then the equation of state is P ~ and the 
differential equation for the perturbations during this period is 

6 



T 



c 2 s k 2 - — U fe = , (72) 



where the conformal time during the matter domination period is r = 2/T-L, corresponding to a cx r 2 , and the velocity 
of sound c s is practically zero. In the limit, c s k — ¥ 0, the two independent solutions for Uk behave like a power law, 
oc t 3 and t -2 . Then the solution for VE^ (oc p l ^ 2 Uk) is given by 

(r) ~ C k + D k r- 5 , (73) 
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FIG. 1: The figure on the left shows V e ff through a transition into radiation domination as a function of r in H e 1 units for two different values 
of q in power law inflation. The transition time corresponds to 5 — 0.01. The figure on the right shows the evolution of V 2fc 3 | ^fc| in A.iz units for 
4 different modes with a = k/a e H e — 5, 1, 0.1, 0.01, and the higher frequency modes corresponding to larger wavenumbers k. The parameters are 
the same as the left figure. 




FIG. 2: On the left hand side, we show the scale dependence of the amplitude of ^ m in H e /Mp units, defined as y/2k^\^ m (k) | , for a transition 
into radiation domination. The two curves are evaluated for S — 0.01. The comoving wavenumber k is in H^ 1 units. The dashed line represents the 
maximum of the analytic expression in Eq. Il67t ■ On the right pannel, we show the value of the maximum for the resonance, ^res, as a function 
of q, for S — 0.05 , 0.01 and 0.001. A larger 5 results in a smaller value of ^res- The solid lines represent the interpolation of the numerical data. 
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where the integration constants depend on the comoving scale k. Thus as the universe enters a stage of matter 
domination, perturbations on scales for which c s kr is small reach certain constant values a few Hubble times after 
the end of inflation. 

As in the preceding section, we consider that the equation of state changes very sharply from the inflationary value 
at the end of inflation, determined by the slow-roll parameter e# in the case of power law inflation, to the matter 
domination value P = 0. The behaviour of the perturbed quantities \& and 1Z during this transition is similar to the 
case of the transition into radiation domination for k <C Ti. c /b. On these scales, we can neglect the gradient terms in 
Eq. (|25l) and therefore the perturbed quantities ^ and 1Z do not grow appreciably during the transition if 5 < en- 
However in the opposite regime, k 3> rl c /S, if we assume that the velocity of sound c s (r) drops down to a negligible 
value during the transition, the potential \ff oscillates with a decreasing frequency, given by W(k, r) in Eq. until 
W(k, t) reaches zero to eventually become negative. At this point, *]/ starts to grow and eventually settle down to a 
constant value during matter domination. 

Assuming that \1/ and 1Z do not vary during the transition, which is valid for k <C Jie/S and 5 < ejj, and taking 
the slow-roll limit \P <g; 72. during inflation, we can determine the integrating constants Ck and Dk- Setting the scale 
factor at the end of inflation to a c ~ 2/r e H e , the potential 4" is given by 



* fc (r) ~ ~H k (r e ) 



1 



(74) 



Thus, a few Hubble times after the end of inflation "J reaches the constant value ^k — 37?.fe(r )/5 on these scales. This 
is the standard result for superhorizon scales during matter domination. Here we have shown that this result will hold 
for k -C H e /5 as long as S < e#. For perturbations on scales k > rl c /S, VP and 1Z vary during the transition because 
the wavelength of the perturbations is the order or smaller than the transition time, and therefore the estimates above 
are no longer valid. Nevertheless, ^ also reaches a constant value during matter domination which in this case is scale 
dependent. 

Approximate analytical solutions are harder to obtain than in the case of a transition into radiation domination for 
k 3> 6 /H c because the gradient term vanishes during matter domination. Neglecting the derivatives of the equation 
of state, the evolution equation for Uk can be written as 



u k = 0, (75) 



where both c 2 {t) and (1 + w) behave like step like functions with the same width determined by the parameter S. 
We have tried 1st order WKB solutions by using Airy functions to match the solutions between inflation and matter 
domination. Despite this solution gives the same scale dependence as our numerical results, that is ^k ~ (k/ric) 1 ^ 2 , 
the amplitude is about a factor of 4 smaller. This is probably due to the fact that the WKB approximation fails 
marginally for large r if c 2 — 0. An asymptotic matching using higher order WKB solutions would reduce the error. In 
Appendix^ we show the details of the first order WKB approximation. The solution with the appropriate boundary 
conditions has the form, 

V6 1/2 f H e 



* fc (r,)~ 4.4 (r^ c r V a/ [jfj • ( 76 ) 

where a = k/He, and t* is the turning point of the differential equation. The turning point marks the time at which 
the frequency of the differential equation changes sign. In this case it is given by W(k, r*) = with 

W 2 {k,r) = k 2 c 2 s {T)-^H 2 (l + w) . (77) 

Numerical solutions of the differential equation in Eq. (f^5|) are shown in Figs. [3] and HI We have modeled the 
time variation of P — P{p) and c s (r) as in the previous section by an error function with a transition time interval 
of At ~ 5/% until they settle down to the values corresponding to matter domination. The numerical results for 
the perturbed quantity are in good agreement with our expectations based on the qualitative analysis. In the 
paragraphs below we describe our numerical results summarized in Fig. [3] and 21 

For large values of q, the superhorizon modes do not show substantial scale dependence and ^ grows and converge 
to an approximate value of ^ ~ 3A-r/5- As q decreases the scale dependence of these modes gets larger increasing 
the red tilt of the spectrum, as shown in Fig. [31 The modes corresponding to scales in the range T~L C < k H c /6, 
reach the approximate value, 

(78) 
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taking the asymptotic limit kr c ^> 1 oilZk in Eq. (I5TH) . On these range of scales "ffc increases as the slow-roll parameter 
€h = l/q gets smaller. (See the right panels of Fig. [3] and Fig. |U) 

On the opposite regime, k ^> l-L e /5, the amplitude of ^k grows with scale as a 1 ! 2 . As k decreases there are a series 
of small oscillations in the spectrum followed by a resonance approximately located at k ~ 5. The resonance is 
as in the radiation domination case the result of the effective potential, V^s, on the modes with wavenumbers about 
k ~ H c / 6. The maximum amplitude of the resonance depends on the height and width of the effective potential V c s . 
These are determined by the parameters q (or en) and S. For a fixed value of 6, the potential barrier described by 
V e g defined in Eq. (|^(7)) behaves as V e g oc 1/(1 + w) 2 , increasing its height as e# becomes smaller. As the height of the 
potential barrier increases, the effect on the a fixed mode k of the field u that propagates across it, is to increase its 
amplitude (as it can be observed in the left panel of Fig. We have described in more detail the effect of a simple 
squared potential on the amplitude of ^k in Appendix [Bj Despite its simplicity, the model reproduces qualitatively 
the features observed in the spectrum of ^k shown in Fig. @] 

Numerical results for the maximum of the resonance in the spectrum of ^k are shown in the right panel of Fig f5] 
A good fitting curve in the regime shown in the figure is 

*res (q,S) = A2—[-^)= —— -ji , (79) 



5 \M P/ u xiPc/ 

where A = 0.195. The dependence of the amplitude of the resonance with q is the same as in the radiation domination 
case. However, in this case ^res diverges in the limit of 6 going to zero. This result is not surprising because if we 
consider an instantaneous transition and use the junction conditions (HJ and ([2"|l. then ^k is given by Eq. (1781) and it 
would diverge in the ultraviolet limit. 

In reality, the sound velocity will not be exactly zero in the matter-dominated stage. In the case of a free massive 
scalar field, the perturbation will not grow on very small scales where the wavenumber exceeds the geometrical mean 
of the Hubble parameter and the mass, k/a > (k/a) c ~ \J Hm |24|. Furthermore, if there is a X(f> 4 self-interaction, 
the critical scale is modified drastically to (k/a) c ~ X 1 ' 2 m 2 /Mp [24| . In any case, these scales provide a natural 
ultraviolet cut-off and regulate the ultraviolet divergence. 

To conclude this section we would like to remark the similar and different features of the solutions in this case com- 
pared to the transition into radiation domination we studied in the previous section. As in the radiation domination 
case, the potential "J does not depend on the slow-roll parameters in the limit k 3> H c /S. It is only in the opposite 
regime, when 7Z does not vary considerably during the transition, that its value at horizon crossing at the end of 
inflation is relevant for the perturbed quantity "J. The scale dependence of the solutions for k < S/H c is similar as 
well, with large values of q giving a nearly scale invariant spectrum while increasing the amplitude of the resonance. 
In contrast, although the height of the resonance is proportional to An, it depends strongly on 8. 



VI. DISCUSSION AND CONCLUSIONS 



We have studied the behaviour of the primordial perturbations through a sharp transition from inflation into the 
radiation and matter dominated epochs. We have found that for transitions that occur much faster than a Hubble 
time, there is a range of scales inside the horizon for which the amplitude of the primordial perturbations is enhanced, 
relative to the amplitude of perturbations that exit the horizon a few efolds before the end of inflation, if slow-roll 
conditions hold during the final stages of inflation. For a Gaussian distribution of the perturbations, this relatively 
large amplitude increase the probability of strong perturbations on these small scales, that may lead to a significant 
production of PBHs. 

For a transition from inflation into a radiation dominated universe, the maximum amplitude of the peculiar grav- 
itational potential is about ~ 2^/3H 2 /tp c . Assuming a Gaussian distribution for the primordial perturbations, 
one can estimate the PBH abundance produced at the end of inflation, using for example Press-Schechter approach. 
(See however [25[ for the effect of non-Gaussian perturbations.) Given that the observed amplitude of the curvature 
perturbation on cosmological scales is roughly given by 1Z ~ 10~ 5 , it turns out that unless the spectrum of TZ is con- 
siderably larger at the end of inflation, the amount of PBH production would not have any cosmological significance 
However, a much larger amplitude on scales k ~ 7i c , is still compatible with current observational data even 
within the slow-roll paradigm [12| ■ Futhermore, PBHs may be overproduced as well with a more complicated scale 
dependence if slow-roll conditions do not hold throughout the whole period of inflation. 

In the case of a transition into a matter dominated universe the production of PBHs could be more dramatic even 
with a very flat spectrum of 1Z on all scales. The maximum amplitude for now depends on delta as ^ oc 1/(5. 
Therefore if the velocity of sound does vanish on a very small scale, such that c s = for k ~ H c /5, there may be a 
significant production of small mass PBHs even if the transition time is not very rapid. 
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FIG. 3: The figures show the time evolution of V 2k s | "l/fc | in A.-jz units through a transition into matter domination for 3 different modes. On 
the left graph, we show the supcrhorizon modes a = k/H c — 0.2, 0.1 and 0.01 for perturbations generated during power law inflation. The smaller 
wavenumbers result in larger amplitudes as a consequence of the red tilt of the spectrum. On the right graph, the corresponding modes arc a — 1, 5 
and 10. In this case the larger wavenumbers result in larger amplitudes. The transition time corresponds to S — 0.01. 
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FIG. 4: The figure on the left shows the matter domination spectrum of ^ , y/ 2k 3 \ | , in H G /M P units for two different values of q. The transition 
time is fixed at 5 — 0.01. The dashed lines represent our analytical approximations onto the left and right of the resonance, 'J'*, oc a. and y/a, 
respectively. On the right, the peak value at the resonance, ^res, is plotted as a function of q for S — 0.05, 0.01 and 0.001. The smaller values of 
S result in larger amplitudes for the resonance. 
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Finally we would like to stress that a detailed analysis of specific models, in which the transition from inflation 
occurs very rapidly, may reveal the overproduction of PBHs, placing constraints on the spectrum of the primordial 
perturbations on scales that are too small for conventional observations. 
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Appendix A: Matter perturbations on scales k 2> He/ '6. 

In this Appendix, we find an approximate analytical solution for the potential \1/ during the matter domination 
epoch on scales such that k H c /5. As explained in Sec.(|V]), we assume that the velocity of sound c s drops down 
from a value of 1 during inflation to zero at the end of inflation much faster than a Hubble time. In this case the 
gradient terms vanish and the simple oscillatory WKB solution we found for the radiation domination case is not valid. 
However, we can still find an approximate asymptotic solution to first order in the WKB expansion that reproduces 
the scale dependence of the spectrum of accurately. 

Neglecting the time derivatives in the equation of state, the evolution equation for the field Uk is given by 



k 2 c 2 (r)~-H 2 (l + w(r)) 



u k = 0, 



(Al) 



where Uk = 2k(p + p) 1 / 2i &k- It is convenient, to introduce an adimensional time parameter r\ = r7i e , where 
H e = H(T e ). Then the equation of motion may be written as 



u = Q{rj)u. 



(A2) 



where Q{rj) is defined as 



' (77) = -% (1 + w)- a 2 c 2 e (77) 



(A3) 



and a — k/a e H e and % = H/H e - 

This equation has a turning point Q (77* ) = at 77* . For 77 < 77* , Q < and the solutions are oscillatory. For 77 > 77* , 
Q > and the two independent solutions may be given by growing and decaying modes. Then during inflation and 
matter domination, the 1st order WKB solutions are given by 



ut f (v) ~ |Q(77)|- 1/4 exp 
uT'iri) - g(77)^ 1/4 exp 



±7 



± 



V\QW\dv 



(A4) 
(A5) 



To match asymptotically the solutions above, we can expand Q(rj) to first order about the turning point and 
use the large argument expansions for the Airy functions. However, it is simpler to apply the following Liouville 
transformation, 



1/2 



Q(v) 



(A6) 
(A7) 



where £ and 77 are analytic functions of each other at the turning point. Then we get the differential equation for Y, 

d 2 Y 



d£ 2 



(A8) 
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where the function # (£) above is 



dt] 



(A9) 



The differential equation for Y becomes an Airy differential equation if we neglect An approximate solution for 
u then is given by 



1/4 



[dAi (0 + C 2 Bi (0} 



(A10) 



which we use below to match the first order solutions across the turning point. 
During inflation £ < 0, and then we have 



e /2 ^ I "ac s (r?)d for, 

Jr), 



(All) 



where r <C and we have taken c s (r) = 1 during inflation. In the models we consider here the velocity of sound 
during inflation always has that value, so from now on we just consider this case. In the limit r) <C r\*, the asymptotic 
expressions for the Airy functions in this limit are 



Ai(-z) ~ n-Wz-V'sm V/ 2 + tt/4) , 



Bi(-z) ~ 7r- 1 /^-V4 C0S g^3/2 + 7r/4 ^ _ 



Taking the initial Bunch-Davies vacuum value for the field Uk, 



2k 



1/2 



-ikr 



one gets the following values for the integrating constants, 



Ci 
C 2 



1 

2k 
-Ci. 



1/2 



t/4 



(A12) 
(A13) 

(A14) 



(A15) 
(A16) 



On the other hand, during the matter domination period £ > 0, and the asymptotic behaviour of the Airy functions 
in this regime is 



Ai(z) ~ 2- 1 7T- 1 /2 z -l/4 e -i 



Si (z) 



2,3/2 



-l/2,-l/4 e §, 3 / 2 



Taking c s (r) and wfr) to be zero for r\ ^> 77*, Q ' is approximately given by 

1/4 



(A17) 
(A18) 



(A19) 



and therefore the solution for VE' during the matter domination epoch, neglecting the decaying mode, is approximately 
given by 



(fi + v) 1 ' 2 



4.4 
TIP 



where we have taken t c = 2/"H e and 



Ml 



~ 8^ 



■He 



(A20) 



(A21) 
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During the matter dominated epoch, the exponential term in Eq. (IA20[) is 



*3/2 



" /3 



-H 2 {1 + w(t)) - a 2 c 2 (r) dr] . 



By integrating by parts expression Eq. (IA22[) , we get 



^ 3/2 = V6 



In (v) ~ \n(r])g- 1/2 —dri 



(A22) 



(A23) 



where q is 



'(v) 



{l + w) 



2k 2 cj 

m 2 



Neglecting the variation of the logarithmic term in the integrand , £ 3 ' 2 is approximately given by 



-e' 2 ~ V61n 
3 \rj 



and therefore the time variation of ^ is negligible, 



^ffc - 4.43 x a 1/2 



He 



7? -0.05 r? -V6 ) 



if r\ is not too large. 

This approximate solution is the first order solution in the WKB expansion. 



(A24) 



(A25) 



(A26) 



The scale dependence of on this 



regime is a 1 / 2 = (k/Hc) 1 ^ 2 , the same that we have found on our numerical estimations. The solution is not time 
independent, although the time dependence is very weak. On the other hand, the approximate value of the amplitude 
in Eq. (|A26I) results in a smaller value than the one we have obtained numerically. These two differences in the 
analytic result, suggest that a higher order WKB solution would correct this deficiency. In fact, for the 1st order 
WKB solution to be valid uniformly for all r it is necessary for Q(r) to decrease much more rapidly than 1/t 2 when r 
goes to infinity [26j . During matter domination Q ~ 1/t 2 , and therefore we believe that a second or higher order WKB 
asymptotic matching would yield a more accurate result. Here we are mostly interested in reproducing analytically 
the scale dependence of the spectrum, and when we compare the analytic with the numeric results in Fig. @J we just 
regulate the lower limit of the integrand in Eq. (|A22[) to obtain a reasonable estimate. We believe that this approach 
is justified given that a relative small error in £ could in principle result in a significant error in * oc exp(2£ 3 / 2 /3). 



Appendix B: The effect of a potential barrier on the potential the broad resonance. 

In this Appendix, we explain qualitatively the origin of the broad resonance and small oscillations that we have 
obtained in our numerical calculations for the spectrum of \I f . In the two transitions from inflation studied here, 
the broad resonance is approximately located at the comoving scale k ~ He/ft and it is followed by a series of small 
oscillations with an amplitude that is strongly suppressed as k increases (see the left panels of Figs H] and UJ). About 
those scales, the derivatives of the equation of state are not negligible, and the function q 2 (k,r) has a complicated 
form that makes difficult to find approximate analytical solutions. For a step like equation of state with a width S/H c , 
these terms typically result in a sharply peaked function with an approximate amplitude of H 2 /5 2 . This suggests that, 
on scales about k ~ H e /S, the terms proportional to w' in q 2 {k, t) act as a potential barrier for the wave equation Uk 
that propagates from the inflation to the radiation and matter domination epochs. 

In order to reproduce the broad resonance in the spectrum then, we adopt a square potential acting as a potential 
barrier, for which solutions can be easily found. The model is very crude, but as we will see it captures the fun- 
damental phenomenon occurring on such scales. We restrict this analysis to the transition from inflation to matter 
domination. This case is easier to analyze because the modes of the potential reach a constant value. The radiation 
domination modes after the transition from inflation oscillate with a decaying amplitude. The quantity of interest 
then is the maximum amplitude during the osillations, which is more cumbersome to calculate. Ncvcrthless, the 
relative enhancement of the modes is caused essentially by the same phenomenon as in the matter domination case. 

We are interested in the evolution of the field Uk after the transition into radiation or matter domination on 
scales such that k > H c /8. We assume that the velocity of sound, remains constant while the wave crosses through 



17 



the potential barrier. This is a good approximation for the radiation domination case because c s does not vary 
considerably during the transition. In the matter domination case, we can consider the idealized situation in which 
c s drops down to zero once the wave function has crossed the potential barrier. To find an approximation solution 
for the perturbations during this epoch then, we can take the resulting solution as the initial condition for (Eq. IA1I) . 

With the simplifying assumptions considered above in mind, we now study the behaviour of the fluctuations in the 
field life for a square potential. The evolution of the field then is determined by the differential equation, 

u{J + [a 2 - V (t)] u k = , (Bl) 

where a = k/H e , and V(r) is given by 

{0, t < -a, 
v 2 , -a<r<0, (B2) 
0, r > 0. 

Taking the initial value of the field u k corresponding to the Minkowski vacuum, 

u k = -^=e~ ikT , (B3) 



/2k 

the solution for r > is 

1 



u k = -±= (Ae~ ikT + Be lkT ) , (B4) 
2k 



where q = \Jk 2 — v 2 , and the constants A and B are given by 

A = {k + qf -(k-qfe 



Akqe- %i ~i- k ) a 
v 2 (l~ e - 2lqa ) 

B = -77 , m ■ (B5) 

As in the previous section, an approximate asymptotic solution valid for throughout the transition and the matter 
domation epoch is 

t \ 1/4 

s ' [D x Ai{i)+D 2 Bi{i)\ , (B6) 

where D\ and D2 are integrating constants. Taking the asymptotic expressions of the Airy functions for r -C in 
Eqs. (|A17I) and (|A18|) . we find that the initial conditions in Eq. (|B4|) are satisfied if 

Dl ~ -^- 1/2 « 1/2 (B-A), (B7) 



D 2 ^Tr^a 1 / 2 (A + B) . (B8) 



The asymptotic solution for r ^> during matter domination for the potential is approximately given by 

where a — k/H c , 77* = r*i? e denotes the time at the turning point (see Appendix |A"|) . and F = A + B. The function 
|-F| which modulates the solution ^ cx a 1 / 2 we have found in Appendix [Al can be written as 



, ,„ a 2 — v 2 cos 2 {Wa 2 — v 2 ) 

\F\ 2 = (BIO) 



for a > v, which is the region of the comoving scale that we are intersted in reproducing. 
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FIG. 5: The fi gurc on the loft wc compare the matter domination spectrum of »3/fc we have obtained numerically with the anlytical approximation 
in Appcndiccsf^and^] The numerical solution corresponds to 5 — 0.01 and q — 100. On the right of the resonance, wc have found a reasonable 
fit with the anlytical approximation for v — 118, I — 0.025 and — 0.2. The dashed line represents our analytical approximation onto the left of 
the resonance, \I/ oc a. 



In the limit k ^> H c we get the trivial result |F| ~ 1, which corresponds to A ~ 1 and B ~ 0. On this regime, 
the variation in the potential does not affect the propagation of the field Uk, and 'Ffc oc a 1 / 2 . As the comoving scale 
k decreases \I/fc decreases until k approaches the scale close to the height of the potential a = v. About that scale, 
the amplitude of "Ffc increases achieving its maximum amplitude at scale about a ~ v. For a < v, the amplitude of 
the perturbations decrease again. The overall effect then is that we find a broad resonance on the spectrum of ^>k 
approximately located at the scale a ~ v. 

In Fig. [5l we compare the result of our numerical calculation of the spetrum of 'Ffc with the analytical approximation 
above in Eq. (|B9[1 . The numerical result is for S = 0.01 and q — 100. To the left of the resonace, we have plotted 
the analytical approximation for k < H c in Eq. (|78p which is a good approximation on that scales. To the right of 
the resonance, on scales such that k > H c /S, we show the analytical approximation above in Eq. (|B9|) for v = 118 
and 5 = 0.02. To regulate the integral in Eas. (|A22l) and (IA23[) we have used the lower limit of integration 77* = 0.2, 
which gives the right result for the scales shown in Fig. [SJ For these values of the parameters v and 5, we see that the 
square potential reproduces the resonance we have found in our numerical calculations sufficiently accurately given 
the crude analytical approximation we have done. 



[1] W. Israel, Nuovo Cim. B 44S10, 1 (1966) [Erratum-ibid. B 48, 463 (1967 NUCIA,B44,1.1966)]. 

[2] N. Deruelle and V. F. Mukhanov, Phys. Rev. D 52, 5549 (1995) [arXiv:gr-qc/9503050l . 

[3] J. Martin and D. J. Schwarz, Phys. Rev. D 57, 3302 (1998) |arXiv:gr-qc/9704049|. 

[4] D. H. Lyth, K. A. Malik, M. Sasaki and I. Zaballa, JCAP 0601, 011 (2006) arXiv:astro-ph/0510647". 

[5] I. Zaballa, A. M. Green, K. A. Malik and M. Sasaki, JCAP 0703, 010 (2007) arXiv:astro-ph/0612379 . 

[6] B. J. Carr and J. E. Lidsey, Phys. Rev. D 48, 543 (1993) [arXiv:astro -ph /9405 027l~ 

[7] B. J. Carr, J. H. Gilbert and J. E. Lidsey, Phys. Rev. D 50, 4853 (1994) |arXiv:astro-ph/9405027] . 

[8] H. I. Kim and C. H. Lee, Phys. Rev. D 54, 6001 (1996). 

[9] A. M. Green and A. R. Liddle, Phys. Rev. D 56, 6166 (1997) [arXiv:astro-ph /970425I1 



[10] S. M. Leach, I. J. Grivell and A. R. Liddle, Phys. Rev. D 62, 043516 (2000) arXiv:astro-ph/0004296 



[11] E. Bugaev and P. Klimai, arXiv:astro-ph/0612659 
[12] K. Kohri, D. H. Lyth and A. Melchiorri, JCAP 0804, 038 (2008) arXiv:0711.5006 [hep-ph]]. 
[13] H. V. Peiris and R. Easther, JCAP 0807, 024 (2008) arXiv:0805.2154 [astro-ph]]. 
[14] E. Bugaev and P. Klimai, Phys. Rev. D 79, 103511 (2009) [arXiv:0812.4247l [astro-ph]]. 

[15] A. S. Josan, A. M. Green and K. A. Malik, Phys. Rev. D 79, 103520 (2009) arXiv:0903.3184 [astro- ph.CO]]. 



19 

[16] L. Alabidi and K. Kohri, Phys. Rev. D 80, 063511 (2009) |arXiv:0906. 13981 [astro-ph.CO]]. 

[17] J. M. Bardeen, Phys. Rev. D 22, 1882 (1980). 

[18] H. Kodama and M. Sasaki, Prog. Theor. Phys. Suppl. 78, 1 (1984). 

[19] M. Sasaki, Prog. Theor. Phys. 70, 394 (1983). 

[20] D. H. Lyth and E. D. Stewart, Phys. Lett. B 274, 168 (1992) 

[21] E. D. Stewart and D. H. Lyth, Phys. Lett. B 302, 171 (1993) [arXiv:gr-q c/9302019 . 

[22] M. Sasaki, Prog. Theor. Phys. 76, 1036 (1986). 

[23] A. R. Liddle and D. H. Lyth, Cambridge, UK: Univ. Pr. (2000) 400 p 

[24] Y. Nambu an d M. Sasaki, Phys . Rev. D 42, 3918 (1990). 

[25] J. C. Hidalgo. larXiv:0708.3875l [astro-ph], 

[26] CM. Bender and S.A. Orzsag, " Advanced mathematical methods for Scientists and Engineers. Asymptotic methods and 
perturbation theory", Spinger 1978. 



